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We present results of non-local and three terminal (3T) spin precession measurements on spin 
injection devices fabricated on epitaxial graphene on SiC. The measurements were performed before 
and after an annealing step at 150 °C for 15 minutes in vacuum. The values of spin relaxation length 
L s and spin relaxation time r s obtained after annealing are reduced by a factor 2 and 4, respectively, 
compared to those before annealing. An apparent discrepancy between spin diffusion constant D s 
and charge diffusion constant D c can be resolved by investigating the temperature dependence of 
the ^-factor, which is consistent with a model for paramagnetic magnetic moments. 

PACS numbers: 72.25.-b,72.80.Vp,85.75.-d 



Apart from its prospects for electronic devices^ sin- 
gle layer graphene (SLG) is also a very promising can- 
didate in the field of spintronics because it is expected 
that spin information can be passed in graphene over 
long distances^ due to the weak spin-orbit coupling and 
low hyperfine interaction^. Up to now, however, the 
measured spin lifetimes in exfoliated SLG (0.5 ns at 
RT— , & 1 ns at 4 K— ) and also in bilayer graphene 
(« 2 ns at RT 7 ) on Si02 are still one order of mag- 
nitude smaller than in conventional semiconductor het- 
erostructures. Even if the mobility fi for graphene on 
Si02 is modified by e.g. ligand-bound nanoparticles 8 
(2700 cm 2 /Vs - 12 000 cm 2 /Vs) or by using high qual- 
ity suspended graphene devices^ (ja > 100 000 cm 2 /Vs), 
measured spin lifetimes are below 2 ns. Similar values of 
r s , slightly over 2 ns, were also reported for graphene epi- 
taxially grown on a semi-insulating silicon carbide (SiC) 
substrat o 10 ! 11 using a direct non-local measurement^ 
while a huge r s was obtained by an indirect 13 - method. 
In Ref. [T2] also a large difference between D c and D s 
was reported, which was first ascribed to spin storage in 
localized states and later reinterpreted using a modified 
^-factor—. McCreary et ali^ studied the influence of 
artificially created paramagnetic moments on spin trans- 
port in graphene, and introduced an effective exchange 
field model leading to an enhanced ^-factor. This vari- 
ety of different results and ambiguity in their interpreta- 
tion motivates further experiments on epitaxial graphene 
to understand the spin relaxation mechanism in order to 
control the spin information for future spintronic devices. 

Here we also use epitaxial graphene grown on the Si 
face of SiC and present non-local and three terminal^ 
spin precession measurements. The latter probes the spin 
accumulation^ 7 - directly underneath the injector elec- 
trode induced electrically by a spin polarized current. 
We compare the results before and after an annealing 
step and observe that our measurements after annealing 
can be well explained with an enhanced ^-factor assum- 
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FIG. 1. (Color online) (a) Schematic drawing of the non- 
local and 3T measurement setup which is used for electrical 
spin injection and detection. The magnetic field B 2 is applied 
along the z axis, (b) SEM picture of the epitaxial graphene 
spin valve device. 



ing that charge (D c ) and spin diffusion constants (D s ) are 
equal. As the temperature dependence of this increased 
factor shows a clear 1/T (paramagnetic) behavior, we 
believe that annealing creates local magnetic moments 
via e.g. vacancies 18 which influence the spin transport 
properties. 

Fig. QJa) shows a sketch of the applied measurement 
methods, Fig. [U^b) a SEM picture of the used epitaxial 
graphene spin valve device. The graphene stripes having 
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FIG. 2. (Color online) (a) Hanle spin precession measure- 
ment (background removed) with a DC current of +10 fx A 
at 1.7 K and fit (continuous line) in non-local configuration, 
(b) 3T measurement (background removed) and Lorentzian 
fit (continuous line). Both measurements are done before an 
annealing step. 



a width of W = 30 /im and a length of about 750 /im 
are produced using a negative resist based electron beam 
lithography (EBL) step and oxygen plasma etching for 
30 s (30 mTorr 2 , 50 W). Afterwards a thin tunnel- 
ing barrier (A\O x ) with a thickness of about 1 nm was 
produced by depositing Al atoms over the entire cooled 
sample (180 K) in a UHV system (jp « 10 -9 mbar) and 
subsequent oxidation in the load lock in pure oxygen at- 
mosphere (p « 3 x 10 -2 mbar) at RT for 30 minutes. 
This AlO^ tunneling barrier with a contact resistance 
R c > 2 kQ provides high spin injection efficiencies and 
reduces spin relaxation induced by the contacts^. The 
ferromagnetic (FM) cobalt electrodes (Co 20 nm) with a 
width of 200 nm (contact A) and 500 nm (contact B) and 
the non-magnetic palladium contacts (Pd 80 nm) were 
each patterned using a positive PMMA resist based EBL 
step. The evaporation is done via electron gun (Co) and 
thermally (Pd) at a base pressure of about 5 x 10 -7 mbar 
followed by a standard lift-off technique. The distance L 
between the edges of the FM stripes is 2 /im. Finally 
the sample is glued into a chip carrier and the measure- 
ments are done using a standard DC setup in a Cryogen- 
ics He-4 cryostat (T = 1.6 .. . 300 K) equipped with a vec- 
tor magnet {B XyVyZ = —1 ... 1 T). The complete sample 
fabrication is done without applying a high temperature 
cleaning step. 

In Fig. [2] typical Hanle curves in the non-local and three 
terminal setup at 1.7 K are shown when using contact A 
as an injector. The FM stripes are magnetized in parallel 
configuration and the magnetic field B z is applied out-of 



plane which leads to dephasing of the spin signal. The 
continuous curve for the non-local curve in Fig.[2fa) is the 
numerical fit with the solution of the following equation 20 
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where P is assumed to be the same for both FM stripes 
and xi, X2 are the points of injection and detection, re- 
spectively. P is the spin injection efficiency, / the in- 
jection current, ujl = g^BB z /h the Larmor frequency 
with the Lande factor g, R s is the sheet resistance of 
graphene, W the width of the graphene stripe, D s is the 
spin diffusion constant, finally r s and L s = \/D s r s are 
the spin relaxation time and length, respectively. The 
Hanle signal Rst in 3T configuration (Fig. E^b)) decays 
with B z due to spin dephasing with a Lorentzian shaped 
according to: 
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As the amplitude of the non-local signal is determined 
by the product of P 2 and L s these parameters are not 
independent in the fitting procedure. Therefore L s = 
1.2 /im is estimated assuming an exponential decay of 
the spin signal given by Rst(B z = 0) at L = /im 
and R n i(B z = 0) at L = 2 /im. Now P and r s are 
the free fitting parameters and D s can be calculated via 
D s = L 2 Jt s . 

In Fig. [2] one can see that for both non-local and 3T 
spin precession measurements the results for P, r s and D s 
are almost identical. This agreement indicates that these 
signals originate from an induced spin accumulation into 
graphene. Slight differences especially in the spin injec- 
tion efficiency P can be explained by a small anisotropic 
magneto resistance contribution of about 0.5 Q of the 
FM stripes to Rst determined in reference measurements 
(not shown) . This small deviation leads to an absolute er- 
ror of L s of about AL S = 200 nm. The spin diffusion con- 
stant D s = 171 cm 2 /s has the same order of magnitude as 
the charge diffusion constant D c = 120 cm 2 /s extracted 
from another epitaxial graphene sample grown with iden- 
tical parameters. This reference sample was measured in 
van der Pauw geometry to determine the sheet resistance 
R s = 4.15 kQ. From low field Hall measurements we get 
the charge carrier density n = 5.25 x 10 11 cm -2 and the 
mobility ji = l/(R s ne) = 2870 cm 2 /Vs. The reference 
sample also showed the quantum Hall effect starting at 
approximately 4 T, underlining the high sample quality. 
Due to a different sample preparation (AlOx on top of 
graphene) for our spin valve device (R s = 1.5 kf2) D c can 
differ from the reference sample and we therefore take D c 
to be identical to D s as determined from the Hanle spin 
signals. By fitting the non-local and 3T measurements, 
we obtain r s = 80 ps which is slightly smaller than in 
exfoliated SLG&±2^. 

In order to check if annealing influences the charge 
transport properties and/or the induced spin accumula- 
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FIG. 3. (Color online) (a) Hanle spin precession measure- 
ment (background removed) at 1.7 K and fit (continuous line) 
in non-local configuration, (b) 3T measurement (background 
removed) and Lorentzian fit (continuous line) . Both measure- 
ments are done after an annealing step. The fittings are done 
treating the ^-factor as a free parameter. 



tion, a post annealing step is done at 150 °C for 15 min- 
utes in vacuum to avoid intercalation of hydroge n 22 ! 23 
via forming gas. Then we repeat the same spin preces- 
sion measurements as before and interestingly we observe 
that r s increases whereas D s is decreased by almost a fac- 
tor of 5 if we assume the same ^-factor as before annealing 
(9 — 9o — 2). For the configuration with contact B as 
an injector we observe an even bigger decrease of D s by 
a factor of about 10 (Table I). At this point we conclude 
that annealing affects the spin transport properties and 
we observe the same apparent reduction of D s as in Ref. 
[l2l where a high temperature annealing step was included 
in the sample preparation procedure. The fact that both 
L s and the 3T amplitude (1.3 Q) at zero magnetic field 
decrease by almost the same factor indicates that the ap- 
plied post annealing step also affects the induced spin ac- 
cumulation underneath the injector electrode. A slightly 
higher sheet resistance (R s = 1.7 kfi) indicates that the 
charge transport properties are only weakly influenced 
after annealing, that means we have the following sit- 
uation: D b c e f° re « Df ter > Df ter . L s = 594 nm is 
reduced by a factor of 2 after annealing and is again ex- 
tracted from the exponential decay of the spin signal at 
the injection point (Rst(0)) and at a distance L = 2 /xm 
OMO)). 

In an attempt to understand the discrepancy of D c 
and D s Maassen et al. 12 first considered localized states 
in the electrically inert buffer layer— (BL) which could 
provide hopping sites for the spins being able to change 
the spin transport properties but not the charge trans- 
port properties. The difference in D s and D c was also 
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TABLE I. D s [cm 2 /s], r s [ps] and ^-factor before and after 
annealing for injector contact A and B. After annealing the 
measurements can also be fitted with an enhanced g e ff > go- 



recently discussed by McCreary et al^. They assume a 
formation of local magnetic moments by Ar sputtering or 
from hydrogen adatoms on exfoliated graphene samples 
which provide an enhanced magnetic field for the diffus- 
ing spins, which can be modeled by an effective ^-factor. 
Maassen et al. 14 reinterpreted their experiments 12 using 
a model of localized states, where the effective Larmor 
frequency is increased in the limit of strong coupling, 
which again can be expressed by an enhanced ^-factor. 

For this reason we also fit our non-local and 3T data 
after annealing, treating the ^-factor in the Larmor fre- 
quency (Eqs. (pp) and fl2J) as a free parameter, and as- 
suming Df ter = Df ter = 171 cm 2 /s. Fig. [3] shows 
that our data can also be well fitted with an enhanced 
^-factor of 6.5 in the non-local and 6.7 in the 3T setup, 
respectively. The oscillations observed in the Hanle curve 
(Fig. EJa)) at higher magnetic fields are phase coherent 
contributions and vanish at higher temperatures. If we 
summarize our experimental findings so far (Table I) we 
can conclude that our measurements after annealing can 
be explained either by D^ ter ^ D®f ter or by an effective 
Lande factor g e ff > 2 as both models reproduce the data 
equally well since Eqs. (pQ) and (j2j) contain the ^-factor 
implicitly in the Larmor frequency ujl and are invariant 
under a rescaling of g, r s and D s . 

To determine which model (hopping or magnetic mo- 
ments) is appropriate in our situation we study the tem- 
perature dependence of spin transport. We observe that 
the enhanced effective ^-factor, as well as the amplitude 
of the spin signal, decrease with increasing temperature 
(Fig. |4j). If g e ff originates from magnetic moments then 
its temperature dependence can be described by the fol- 
lowing equation^: 
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This is the low field approximation of the Brillouin func- 
tion of a spin 1/2 paramagnetic material. A ex is the 
strength of the exchange coupling, tjm represents the fill- 
ing density of the magnetic moments, go = 2 is the g- 
factor for free electrons and ks is the Boltzmann con- 
stant. As it is seen in Fig. [5] the measured temperature 
dependence of g e ff is well described by Eq. (j3j). This 
temperature dependence is compatible with the effec- 
tive exchange field model proposed by McCreary et al^ 
which describes the enhancement of the magnetic field 
felt by the diffusing spins due to localized paramagnetic 
moments. Maassen et ali^, on the other hand, interpret 
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FIG. 4. (Color online) Temperature dependent non-local 
Hanle measurements (background removed) with ^-factor as 
further fitting parameter. 



their data by hopping of the diffusing spins into localized 
states, whereby they experience a dramatically higher 
g-factor. In their work the enhancement of g is most 
pronounced at room temperature, whereas in our case 
the maximum g e ff is obtained at low temperature. We 
therefore believe that post annealing creates an amount 
of randomly positioned magnetic moments resulting in 
an increased effective magnetic field B e ff = B z + B ex 
composed of the applied out-of plane magnetic field B z 
and of an exchange field B ex coming from the induced 
magnetic moments. This enhanced magnetic field can 
be modeled by an effective ^-factor in the Larmor fre- 
quency ujl- As we nearly get the same Lande factor 
from the non-local Hanle fit and the Lorentzian fit we 
can not decide if the magnetic moments are formed in the 
graphene/buffer layer transition or at the AlO^/graphene 
interface. The difference between our experiment and the 
work of Maassen et al 12 may be due to different anneal- 
ing conditions. 



FIG. 5. (Color online) Temperature dependence of the g- 
factor. The continuous curve is the fit according to Eq. J3}. 

In conclusion, an electrically induced spin imbalance 
from ferromagnetic Co stripes can be analyzed via spin 
precession measurements in both non-local and three ter- 
minal configuration. By introducing a post annealing 
step, we observe that the spin relaxation length as well 
as the non-local and 3T Hanle amplitude decrease. Fit- 
ting of the non-local and 3T data after annealing show an 
increase of the ^-factor if spin and charge diffusion con- 
stants are assumed to be the same. The origin of the g- 
factor enhancement are local magnetic moments formed 
by annealing. The reduced spin lifetime and length sup- 
port this assumption because local magnetic moments act 
as an additional spin scattering source. Finally, the tem- 
perature dependence shows a clear evidence that para- 
magnetic moments are created as the effective ^-factor 
scales with 1/T with increasing temperature. 
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